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AND GENERAL BODY OF REVOLUTION 

By G. Lotns Smith 


SUMMARY 

The electromagnetic field equations are applied to 
nonferromagnelic conducting slowly spinning thin - 
and th ick-wall open-ended cylinders and th in -wall 
cone frustums , and also to thin-wall tumbling cyl- 
inders , to calculate the induced eddy currents from 
which the resulting torques are det r mined. A 
method is also presented for determining the eddy 
currents , and hence the torque , in a series of cone 
frustums joined end to end . When applied to the 
limiting case , this method leads to the solution jor 
the general body of revolution. 

Figures that show the variation of torque with 
fineness ratio and thickness ratio are presented for 
thin- and thick-wall cylinders . The torque acting 
on a tumbling cylinder was found to be one-half the 
torque acting on a symmetrically spinning cylinder , 
all other factors being equal. Results from this 
analysis are directly applicable to calculations of 
the torque acting on spinning and tumbling satellites . 

INTRODUCTION 

One of the many considerations in dealing with 
satellites is the interaction of the earth’s magnetic 
field with the conducting shell of a spinning satel- 
lite. (See refs. 1 to 4.) The geomagnetic field 
induces eddy currents within the rotating con- 
ducting shell, which in turn interact with the field 
to produce a torque. One component of this 
torque slows the rotation, and another component 
tends to precess the direction of the spin axis. 
(See ref. 4.) 

In studying this interaction analytically, it is 


necessary to apply Maxwell’s electromagnetic field 
equations with the appropriate boundary condi- 
tions to the problem. The solution for a sphere, 
which is applicable to spherical shells such as those 
of the Vanguard series of satellites, has been pre- 
sented in references 4 and 5. The purpose of 
this paper is to present the solutions for some 
other configurations. Two cylinder cases are 
solved: the cylinder spinning about its center 
line and the cylinder spinning about an axis per- 
pendicular to the center line. All rotations can 
be resolved into these two components. The im- 
portance of having both solutions available for 
the general cylindrical satellite is that, although 
it may he injected into orbit spinning about (he 
axis of minimum moment of inertia, internal dis- 
sipation of kinetic energy causes the spin axis 
eventually to shift to the axis of maximum mo- 
ment of inertia. 

Thin cones and frustums of cones are also 
investigated and equations are derived for the 
induced currents within the shell, from which the 
torque follows. A method is set up for extending 
these results to a series of frustums joined together. 
This method leads to a solution for a general body 
of revolution. 

In each ease, equations are derived for the cur- 
rent density throughout the body and for the total 
resultant torque. 

SYMBOLS 

C curve of integration 

c velocity of light 

E electric field intensity vector 
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E (with subscript) component of E 


unit vectors along f , r), and f axes 
vector force per unit volume 
magnetic field intensity vet? tor 
magnitude of H 

unit vectors along X, Y, and Z axes 
current density vector 


e o e 2? e 3 

F 
H 

h 

j, k 

J 

J (with subscript) component of J 


JJ ) 

kni) kmi 


L 

/ 

r 

r 

t 

u 

w 

*,y,z 

Yn( ) 

Z *{ ) 


Pm 

e 

x 

M 


f, v> f 


£T 

r 

4> 

(0 

CO 


Bessel function of first kind, of order n 
eigenvalues in thick-wall cylinder solu- 
tion 

torque vector 
length of cylinder 
function defined by equation (36) 
radius vector 

radius (cylindrical coordinate) 
time 

subscript) unit vector, in direction 
indicated 

transformation variable, defined by 
equation (71) 

distance along coordinate axes X, T, 
and Z 

Bessel function of second kind, of 
order n 

cylinder function of order n, E mn J n ( ) + 
E mn Y n { ) 

quantities defined by equation (40) 
quantities defined by equation (46) 
angle from X-axis in X,Y plane (cy- 
lindrical and spherical coordinate) 
angle between Z-axis and H 
angle defined in figure 10 for (um- 
1 )1 ing-eyl i n dor an alysis 
polar angle in plane, used in cone 
analysis 

coordinate axes used in analysis of tum- 
bling cylinder (see fig, 10) 
distance along cone from vertex to 
point (spherical coordinate) 
electrical corn! uet i vi ty 
thickness of thin wall 
harmonic function; for example, equa- 
tion (13) 

cone half-angle (fig. 11) 
stream function 
spin vector 
spin rate 


Subscripts: 

a, b quantity evaluated at end a or end b 

i inside 

o outside 

av average 


GOVERNING EQUATIONS AND BOUNDARY 
CONDITIONS 


In order to formulate the problem mathemati- 
cally, it is necessary first, to list the equations and 
the boundary conditions to be used. Maxwell's 
equations describe the magnetic and electric fields 
inside and outside the cylinder. (See refs. 4 and 
6.) The analysis is restricted to nonferromagnetic 
metals so that, in Gaussian units, the permittivity 
and permeability are near unity. Also, for spin 
rates reasonable for most satellites, the magnetic 
field due to the induced eddy currents is small com- 
pared with the primary field; therefore, the unper- 
turbed magnetic field can he used in the electric 
field equations. This approximation is justified 
in reference 4. Also, for spin rates of the magni- 
tude applicable to satellites, the charge density 
within a conductor will be negligible. Tiius, as is 
shown in reference 4, the electric field equations 
are 


VXE=~c 


dH 

dt 


( 1 ) 


VE=0 (2) 

In stationary axes, the electric field can be 
written (ref. 4) as: 

E=v^+c" I («Xr)XH (3) 

where v 2 $ — 0, and the potential of the electric 
field, is determined by the boundary conditions. 
The term c -I (wXr)XH is the induced electric 
field. The current follows immediately from 


J— crE (4) 

The force per unit volume is then 

F=c'JXH (5) 

and the torque is calculated by integrating the dif- 
ferential torque 


dL— rXdF (6) 

The only boundary condition is that the compo- 
nent of current (or electric field) normal to the 
surface vanish at the surface. This condition, 
together with the induced field, is sufficient to 
determine 4> and completely. 

Equation (3) is well suited for calculating the 
electric field in symmetrically spinning cylinders. 
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inasmuch as the boundary condition can be readily 
applied in this case to determine <f>; it will be used 
in the study of both symmetrically spinning thin- 
wall and thick-wall cylinders. Plowever, in cal- 
culating the electric field in thin -wall tumbling 
cylinders and spinning cones, equation (3) is not 
so well suited, and it becomes convenient to use 
a stream function to solve equations (1) and (2) 
simultaneously. 

SYMMETRICALLY SPINNING CYLINDER 

For the symmetrically spinning cylinder, an 
approximate solution is first obtained for the 
thin-wall shells. The solution for the thick-wall 
cylinder is then derived, and a comparison is made 
of the two solutions. 

To study the case of a symmetrically spinning 
open-ended cylinder, Cartesian and cylindrical 
coordinate systems are first set up as shown in 
figure I. The Z-axis is set up along the center 
line of the cylinder, and the X-axis is defined in 
such a way that H lies in the X, Z plane and forms 
an angle X with the Z-axis. The quantities H, r, 



and w can then be written as 
H— A(i sin X+k cos X) 

=h(u, sin X cos 0— u<? sin X sin 0- \- u* cos X) (7) 

v=u r r+u z z (8) 

(9) 

Equation (3) thus becomes 

E=V$-\-c- l ho)r(u r cos X— u z sin X cos 0) (10) 

The boundary conditions may then be written as 

E : (r,d,±^)= 0 (11) 

E r (r t , 6, z) = E,{r 0 , d, ~)=0 (12) 

THIN-WALL CASE 

For a t lun-wall open-ended cylindrical shell, 
the radial component of flow will be negligible by 
comparison with the circumferential and longi- 
tudinal components, J e and J Z} respectively, and 
Je and J z will not vary significantly between 
r~Ti and r~r 0 . The problem therefore ran be 
considered to be primarily dependent on 0 and 
z. Since the cylinder can then be cut along an 
element and developed onto a plane, the problem 
may be treated as two-dimensional, with rO as 
the abscissa and ^ as the ordinate, where r is 
taken to be the “average” radius. Only a strip 
of the plane one period in width need be considered. 
A potential field Y<£ is now superimposed on t lie 
field and adjusted to make the longitudinal com- 
ponents of the total field vanish at the boundaries 
of the region corresponding to the open ends of 
the cylinder. The potential then is 

oo 

$=VhVl., sin 7/0+ B n cos nB) 

n =0 

(c n sinli ~ +D n cosh ”') (1 3) 

Equation (10) thus becomes 
E=— n^-^rsin X cos0+u* ^2(A n sin 7/0 

7?=0 

+ 7+ cos 7/ 0) (c n cosh y+Z> n sinli ~ 

+~ X n(A H cos nQ—B n sin nB) ( C n sink — 

T n =0 \ r 

+T) n cosh — ^ (14) 

r * 


Fighrk 1 . — Coordinate systems for symmetrically spinning 
cylinder. 
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The radial components have been dropped in 
equation (14) for this thin-wall case; thus, equa- 
tions (12) are satisfied automatically. In order 
to satisfy the boundary condition (11), it is 
necessary that the first term of equation (14) 

cancel the first summation at This is 

Zi 

easily accomplished by letting all the coefficients 
be zero except B x and C x , Then 


BO l 

0——c- l hur sin A cos 0+--— cos 0 cosh ?r- 

r 2r 


from which 


B X C X 


c 'hxor 2 sin A 

cosh “ 

2 r 


(15) 


Substituting equation (15) into equation (13) 
gives the potential as 


sinh 


<f>(0, z)=c sin A cos 0 
The electric field is then 

E(0, z) =— U £ e~ 1 h( J or sin A cos 0 

—U $c~ l ho>r sin A sin 0 


(16) 


cosh 


2 r 


cosh - 
r 


cosh ; 


sinh - 
r 


cosh — 
2r 


(17) 


The electric field having been determined, the 
current follows immediately by equation (4). The 
torque is then calculated by 


L=c 



rX(JXH )dV 


( 18 ) 


where V is volume. Equations (7), (8), and (17) 
are substituted into the integrand of equation (18). 
The resulting vector expression in terms of u r , u<>, 
and u z is then referred to the X,Y } Z system in 
terms of i, j, and k. The result is then integrated 
over the surface of the cylinder. (Because of the 
thin-wall approximations, the integration with 
respect to r is replaced by simply multiplying by 
the thickness t.) The final result is 

L—-7r(TC~ 2 h 2 io sin A Ylr ^1— y- tanh — ^ 

(i cos A— k sin A) (19) 


2? I. , . 

The factor 1 — r tanh — > which is the torque 
L 2 r 

per unit length normalized with respect to the 
torque per unit length of a similar cylindrical shell 
of infinite length, is plotted as a function of fine- 
ness ratio in figure 2. It is seen from this figure 
that the torque per unit length varies rapidly 
with fineness ratio up to a ratio of approximately 
5, after which the torque per unit length is a weak 
function of fineness ratio. The torque, similarly 
normalized, is shown in figure 3. 

The conventional stream function \p f describing 
the current paths within t lie cylindrical shell, is 
defined by 


J — 17 


dz 


— — J e = — cr Ke 


( 20 ) 


Dimensionally, this definition corresponds to a 



Figure 2. — Normalized torque per unit length as a 
function of fineness ratio for spinning thin-wall cylinder. 



Figure 3. — Normalized torque as a function of fineness 
ratio for spinning thin-wall cylinder. 
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unit thickness. Then by equation (17) 


^(0, s) = — ac-'hwr 2 sin X sin 0 



( 21 ) 


The current paths, or streamlines, are given by 
lines of constant ^ and are shown in figures 4(a), 
4(b), and 4(c) for fineness ratios of 2, 4, and 8, 
respectively. The mapping for will be 

identical to that shown for 0^0 ^w. The values 
of ^ have been normalized by dividing by the 



Figure 4 - Current paths for spinning tliin-wall cylinder. 


maximum value. The direction of the flow de- 
pends, of course, on the relative direction of the 
spin vector and the applied magnetic field vector 
and thus is not indicated. Comparison of these 
figures shows physically why the torque per unit 
length varies as it does with fineness ratio. As 
the fineness ratio is increased, the current paths 
become more nearly straight and parallel, except 
near the ends; thus, in the limiting case of an 
infinite cylinder, the streamlines are parallel. 

In the preceding analysis the problem of the 
torque and eddy currents produced by a conduct- 
ing cylindrical shell spinning in a magnetic field 
lias been studied on t he basis of tliin-wall approxi- 
mations. Exact solutions to the problem for a 
thick-wall cylinder will now be derived because 
of their intrinsic interest and also to substantiate 
the thin -wall treatment and to find its limitations. 

THICK-WALL CASE 

Tn order to determine the electric field within a 
thick-wall cylinder, equations (10) to (12) arc 
again employed. The cylindrical harmonic may 
be written as 

$(r,0,2)=Xl S (Amn shill kmnZ + Rmn 

m n 

cosh Jc mn z)(C mn sin nd \-D mn cos nfi)Z n (k mn r) 

(A n sin n6-\ K, t cos nd) (C a z+D n ) (E n r n 

n 

+ F n r~ n ) +C 0 log. r (22) 

where 

Z n (k m „r) =E mn J„ ( k mn r ) + F mn Y„ ( k mn r ) (23) 




Figure 4. — Continued. 


Figure 4. Concluded. 
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and J n { ) and V H ( ) are Bessel functions of order 
n of the first and second kinds, respectively. For 
simplicity, X will be taken to be *72. From 
physical considerations it is apparent that only 
the component of the magnetic field normal to 
the spin axis is effective in generating a current; 
therefore, this restriction on X will be inconse- 
quential. By using equations (10) and (22), the 
components of the electric field vector within t lie 
conductor can be written as 



=E2 (Am,, si nil k mn z 

m n 



~\'B mn cosh ^-mnw)(Cn« sill 1(9 
+ T)„ ul cos nd)Z'Jk„,„r)+y'.(A„ sin n9 

n 

+ B„ cos 7 iff) (C n z +D n )(E„r n - 1 

— F n r~”~ 1 )n + ~- (24) 


= ;SS A„,„ sinh k mn z 

' m n 

+ B„,n cosh k„ m s) (C mn cos nd 
— D mn sin h9)nZ„(k mu r) 

+~ X A« cos 7)0— B n sin nO) (C„z 

* JT 

+ D„) ( E n r a + Fr~*)n (25) 

E t =^—c~ l 1u>or cos 6 


; I mn (^ 1 m n COSll Ic mn Z 

m n 

+B mn sinii 

kmn ~) (C mn sin nd 

+E mn cos »9)Z n (k mn r) 

+S sin nd+B n cos nd) (E n r" 

n 

+ F n r^ n )~c~ l h^r cos 6 (20) 

The boundary conditions are now applied to t lie 
problem. First the requirement of equation (11) 
that the longitudinal component of the electric 
field vanish at the (aids is applied to equation (26). 
By symmetry B mn =0. Also, only the cos 9 terms 
can have nontrivial coefficients, that is, D mn =^ 
B n =0 (n^l) and C mn = A*— 0. Next, the re- 


quirement of equation (12) that the radial com- 
ponent of the electric field vanish at the inside 
and outside surfaces is applied to equation (24). 
This condition gives C o ~0, E tl ~- F n ~0, and 

7J n {k mn r 0 ) = Z;(t w „r < ) = 0 (27) 

Equations (27) yield the characteristic equation 
which determines the eigenvalues k mn and also 
p 

The method of evaluating them will he 

LLr mn 

discussed presently. By letting D m j— 1 and 
dropping the n subscript, equations (22), (24), 
(25), and (26) are reduced to 

<£(r, 6, 2 )= X Am sinh k m z cos dZ x (k m r) (28) 

m = \ 

/?,= 7 A m sinii k m z cos dZ[(k m r) (29) 

m - 1 

Ee=— XI sinh k m z sin 9Z x (k m r) (30) 

7 / 1=1 


E £ = X A m k m cosli k m z cos OZi{k m r)—c“ l hwr cos 6 

771 = 1 


( 31 ) 


F m 


In order to evaluate k m and j^- } equation (27) 

‘ -*mn 

is expanded by means of equation (23) and then 
Z[(k m r) is replaced by k m Z u {k m r)—^ Z x (k m r). The 
following two equations are thus obtained: 

E m ^k m Jo (k rn r t ) — y J x (k a rt) J 

+ F W Tj (& m rj)J=0 (32) 

Ef n tnJ o(k m 1 o) ~ 1 (krn^o) 

+ F m 0 (k m r 0 ) — y 3 i {k m r a ) J“0 (33) 

y 

Eliminating ’ n between equations (32) and (33) 

m 

gives tlie characteristic equation for the k m values; 
k m *J o (I m? o) ~ J 1 (J^nP o) k tH J 0 (Ji. tt\P i) J l {,k /n T t ) 

*o T_ l 

k/n^ u ilnjo) 3 l{k m /J) I m 3 0 {k m V ,•) * 3 i 

Tq *1 

( 34 ) 
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One method of solving equation (34) is to re- 
write it as 


^ in? (,*} 0 (J\. m To) fj t ( J^mTo, ) k m T j*J q {kmT i) J I (kmT i) 

fcmVo^- o(kmVo) ^ rnTo) k^X fi o 1 ( J^nJ i ) 

(35) 

Now define P{x) by 


P(x)- 




(36) 


~xY 0 (x)—Yi ( x ) 

Thus, the characteristic equation may be written 
as 

P(M = P(k m r 0 ) (37) 


A plot of the variation of P(x) with x which is 
applicable for all cases can be made. To deter- 
mine the eigenvalues for a cylinder with a given 

rilr 0 , P (S') is plotted as a function of ,r. The 

intersections of p(^yd£^ with P{x) satisfy equa- 
tion (37) and therefore give the desired eigenvalues 
k m =^ — Such a plot is shown in figure 5. The 
solid line is for P(x) and the dashed line is 

for P(0.5x) or —=0.5. In this example, the first 
To 

three eigenvalues given by the intersections are 
seen to be k m —\A2, 6.53, and 12.65, where r 0 is 
assumed to be unity. The variations of the first 
few eigenvalues with rjr 0 are shown in figure 6. 

After the eigenvalues have been determined, 
F m /E m is given immediately by 

& = -P(k nl r 0 ) (38) 



Figure 5.— Plot for graphical determination of eigenvalues. 
644936—62 2 



Figure C>. — Variation of eigenvalues with rjr 0 . 

The variation of the first few values of F m /E m 
with Vi!r 0 is shown in figure 7. The singular 
points occur where E m becomes zero while F m 
remains finite. 

The A m values are determined by again apply- 
ing equation (11) to equation (31) and using the 
orthogonality properties of Bessel functions: 


A 


m 


c 1 hwr 0 a m 
k m cosh Y 


(39) 


where, by reference 7, 


f ° r 2 Zi (k m r)(lr 
r„J r[Z, (k M r)f<lr 

k m 

r„ tt UZi (k m r)] 2 —Z 0 {k m r)Z,(k m r ) } |'» 


(40) 


The a m values are thus nondimensional and arc 
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Figure 7. — Variation of F m /E m with rjr 0 for different 
eigenvalues. 

functions only of r f /r 0 for a given m. Tt is noted 
that E m and F m need not be separately determined ; 
only their ratio F m jE m need he found. Substitut- 
ing equation (39) into equations (28) to (31) gives 

$(M, z)=c~ 1 h(i)f e cos 6 ± ^ —^.Z x {k m r) (41) 
m ~ lKn cosh =sL 

E r =c~ l hur 0 cos 6 ±, r — HH Z '^ k "‘ r ) (42) 
m=l m cosh^ 


E e =-c~ l hu’f sin 6 JZ £= — 1 - ^ Z, (k m r) (43) 


<—l b b f 

m = 1 *w/t i 


cosh 


7 ?,— — c 1 hcor 0 cos 0 


V ’sr - v C0sllA m 2 *7 (1 \ 

m=l i a fltt- 


cosh 


(44) 


The electric field is therefore determined within 
the cylinder. The convergence of these summa- 
tions is generally good. 


The procedure for determining the torque is 
essentially the same as outlined before for the 
thin-wall cylindrical shell. The result is 

\j— 7 rcrc~ 2 h 2 sin (^) J 

-s ~^~r tanh (k'„ Jr)\( • cos X-k sin X) 

771 = 1 ff * \ 0/ J 

2 r. (45) 

where lj2r u is the fineness ratio, k’ m is the eigen- 
value noiulimensionalized with respect to r 0 , that 
is, k' m =--r a k m , and /3„, is defined by 

A.=£? J"' E7^k m r)dr=°± [£ Z 2 (*„r)][’ 

(46) 

Note that again, for large values of the fineness 
ratio lj2r 0} the torque per unit length is inde- 
pendent of the fineness ratio. 

The torque per unit length per unit thickness 
as a function of r 2 /r c is shown in figure 8 for 
various fineness ratios. (The product ac~ 2 h 2 o)r 0 A is 
taken to be unity.) Figure 9 is a cross plot 
showing the torque per unit length per unit 
thickness as a function of fineness ratio for various 

rJr 0 values. The — =1.0 curve is identical to 

r d 

that for the thin-wall solution shown in figure 2, 



Figure 8. — Torque per unit length per unit thickness as 
a function of rjr 0 for spinning thick-wall cylinder. 
The product <rc~-h-wr 0 A is taken to be unity. 
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Fir; the 9. — Torque per unit length per unit thickness as 
a function of fineness ratio for spinning thick-wall 
cylinder. The product ar -2 /* 2 cor„ 4 is taken to be unity. 


THIN-WALL TUMBLING CYLINDER 

The coordinate systems used in the analysis of 
the tumbling cylinder are shown in figure 10. 
The £, rj, f coordinate system, with the unit 
vectors e l7 e>, and e 3 , is space fixed at the center 
of the rotating cylinder, with the £-axis parallel 
to the spin axis and the f-axis oriented so that 
the magnetic fieljj vector lies in the £, f plane. 
The X, Y , Z syster^ with unit vectors i, j, and k is 
fixed in the cylinder with the origin at the center 
of the cylinder, the X-axis being alined with the 
£-axis, and is rotated from the space-fixed system 
by the angle fj.. Tn addition, a polar coordinate 
system r,d,z with unit vectors u rj id, and u* is 
fixed in the cylinder, 0 being measured from the 
-Y-axis, as in figure 1 . 

Now, because only the normal component of 
VXJ is effective in a thin-wull conductor and 
radial currents are negligible, equations (1) and 
(4) give for this case 

(VXJ)r =^-^7 = “ffC" 1 /t« sin X cos ju sin 6 

(47) 

The continuity requirement (eq. (2)) reduces to 

< 4s > 



Radial current having been neglected, the bound- 
ary condition is simply 

J,= 0 (s=±0 (49) 


This problem could be solved l) 3 T a scalar potential 
as before. However, to demonstrate an alternate 
approach, the solution will be by means of a 
stream function. With a stream function \p de- 
fined as before (see eq. (20)) 

— (VXJ)r (50) 

so (hat, by equation (47), 

V 2 \f/= ^ 2 2 + aC ^ t0 ° S ^ n ^ C0S 0 (51) 

Solving equation (51) yields 
ip=(Tc~ 1 hwr 3 sin X cos ^ sin 0 

CO 

+’S (-4,, sin 7/0 + B cos nB) 

n= 0 

(c„ sinh n j+D n cosh " ') (52) 

Only the sin 6 component of the harmonic part 
of \p will not vanish when the boundary condition 
is applied; therefore, /1„ = 0 (n-> 1) and B n ~0^ 
By symmetry C n ~ 0. The current density is 
then calculated by equations (50) and (52), 
under the condition that the longitudinal com- 
ponent vanish at the ends (eq. (49)), and A x Di is 
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+©3 cos X cos 2 p) (57) 

Averaging tliis torque around one revolution 
gives simply 

L av = ~ <rc~ 2 h 2 <j} sin XYlr tanli 

e A sin X+e 3 eos X) (58) 

Comparison of equation (58) for a tumbling 
cylinder with equation (19) for a symmetrically 
spinning cylinder shows the two expressions to be 
identical except for a factor of 1/2 in the case of 
the tumbling cylinder because of the sinusoidal 
variation of the current. Also, it is seen from 
equations ( 21 ) and (53) that the streamlines are 
iden tical. 

THIN-WALLED CONES AND CONIC FRUSTUMS 

In studying the magnetic torques on thin-walled 
symmetrically spinning cones and cone frustums, 
coordinate systems arc set up as shown in figure 
11 . Cartesian and polar coordinates are oriented 


as before, the Z-axis being parallel to w, the X, Z 
plane containing H, 0 being measured in t lie 
X , Y plane from the X-axis, and r being measured 
normal to the Z-axis. The origin is placed at the 
vertex of the cone. The cone half-angle is <p. Tn 
addition, p is defined as the distance from the 
vertex to a point on the cone, and unit vectors 
14 and Up are defined normal to p and u 6 and 
parallel to p } respectively. 

As before, only the component of VX J normal 
to the surface is considered. Equation ( 1 ) gives 

(VXJ)*=“< 7 c~ l hu> sin X cos <j> sin 9 (59) 

which is the governing equation for a conical 
surface. 

As with the cylinder, the cone or cone frustum 
is a developable surface and can be rolled out on 
a plane; therefore, the problem becomes a bound- 
ary-value problem in a sector of an annulus or 
of a circle in a plane. The polar coordinates in 
the plane are the radial distance p and the central 
angle v , which is related to the angle 6 by the 
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equation 

v=dsin 0 (60) 

(This equation is readily verified by noting that 
a circumferential distance (p sin <t>)d on tlie cone 
becomes pv in the plane.) 

Tn polar coordinates, the Laplacian of the 
stream function is given by 


dp 2+ pdp + 


5 V 

p 2 d? 


so that equation (59) now becomes, with the aid 
of equation (60), 


1 


dfy 


<) 2 0 d0 ^ Y 

dp 2+ pdp + p 2 sin^ W^° c h<J> Sln X cos * sin 0 

(61) 


As in the case of symmetrically spinning cylinders, 
only the sin 0 term of 0 remains after the boundary 
conditions have been applied. 

Accordingly, the solution of equation (61) is 
assumed to be of the form 


Wp> 0)=sin ef(p) 


(62) 


Substituting this expression into equation (61) and 
rewriting gives 

P 2 ^{+p Jp~~ cs ° 2 (TC~ l ho) sin A cos 0p 2 (63) 

The result is seen to be an equid intension al 
equation. Substituting the solution for f(p) into 
equation (62) gives 

^(p, d) = ac s . u i e (p*-\-Ap csc * 

4 — csc 0 

+ 7?p -csc *) (64) 

where /I and 7? are to be determined. 

The two components of current density are 

d0 T , a c~ l hu sin A cos 0 sin 0 

=Je{p } 0)=— . - 2 ^ 

dp 4— esc 0 

(2 p+A csc c />p csc 7? csc 0p‘ c5c ^" 1 ) (65) 

<rc~ l hw sin A cot 0 cos 0 

4— csc 2 0 

(p+ylp C8C t-'+Bp-*™ +- 1 ) (66) 


In applying these equations to a cone, the 


boundary conditions are that J 9 — 0 at the end 
P—Pa and that the current density remains finite 
throughout, the cone. By the latter condition 
B~0, and the former gives 

Pa+Ap / sc *- I =0 
or 

A= — p a 2 ~~ CBC * (67) 

For a cone frustum, it is required that J p = 0 at 
both ends, p—p a and p=p&. Thus 

Pa +Ap a csc ^- I +7?p a - c8c *- ! = 0 (68a) 

p b +A Pb C9C *-' + Bp b ~™ c = 0 (68b) 

Simultaneous solution of these two equations 
gives 

_ - CSC <£- 1 „ „ A -CSC 0-1 

a Pa Pb Pa Pb f fiQ ^ 

- — p a csc «-l pr csc *-J_ pa -csc ^-Ip^csc *-* ^ 




_ cac 0-1. . csc 0-1 _ 

PaPb Pa Pb 


p." e * _1 pr e,e *“ I - 


A —csc 0 — 1 _ cac 0 — 1 
' Pa Pb 


(69b) 


Substitution of the current expressions (eqs. (65) 
and (66)) into equation (18) gives the torque: 

j Trac~ 2 h 2 u) sin A cos 2 0 sin 0r / p 4 Ap C3C * +2 

4— csc 2 0 \4 csc 0+2 

7? -CSC0 + 2X 

+ ^-- 1 ^ (— i cos X+k sin X) (70) 

Tt is seen that the solution in the form given 
here contains a singularity for csc 0— ±2 or 
<£=-^30°. This is due to the homogeneous part 
of the solution becoming identical with the in- 
homogeneous part. The solution for this case 
may be obtained by introducing a transforma- 
tion variable w defined by 


w=l«g< P J 

whence 

df = df 

p dp dw 
p dp 2 dw 2 dw 

With these expressions, equation (63) gives, for 
0^+30°, 


d 2 f u -u * 
- r +—4/~(rc Via? sin 
dw 



e 


2w 
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The solution of this equation is 

(5 

J=^§ ac~ l hw sin \(we iw +Ae iu! +Be~ 2u ’) 

O 

<tc _1 Ao) sin A(p 2 log t . p+^4p 2 +^p -2 ) (72) 

8 

The stream function is thus given by 

4'(p 1 9) = ^ <TC~^ho3 sin X sin B(p 2 log, p~\~A p"-f~ B p *) 

o 

(73) 

and the current density components are 

Je{p, 0)=^r *c~'hu sin X sin 0[p(l +log, p ) 

-\-2Ap — 2 Bp~~ 3 } (74) 

/g 

J P (p, 0) = — vc~ l hu sin X cos 0(p log,, p 

+.lp+B p - 3 ) (75) 

As before, these equations are applicnl to a cone, 
and the constants are found to be 

log, Pu 

B^O 

Likewise, for a frustum, the constants are 

4 . P * 4 log, Pb—Pa log, p a 

Pa ~~Pb 

P«7 4 pt> 4 (l0g r p a — lOg, p b ) 

** 4 4 

Pa Pb 

The resultant torque is then 

L=-^ TT<TC- 2 k 2 o> sin Xr Q p 3 ^log, p— Ap 2 

— ^ i?p~ 2 J (i cos X— k sin X) (76) 

SERIES OF CONE FRUSTUMS 

If a series of in conic frustums are joined end 
to end, as in figure 12, equations (64), (65), and 
(66) apply in each section, and it remains only to 
determine the A and B for each section in order 
to define fully the current and hence the torque. 
At the joint, the radial component of current 
must be continuous, and the circumferential 
component must also be continuous as otherwise 
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FmvRE 12. — Series of cone frustums. 


a vortex line would be formed, with resulting 
infinite curl of current along the joint. Therefore, 
the conditions for determining the constants are 
that the respective components of current arc 
continuous til the joints, and that at an open end 


the radial component vanishes, or that, if the 
end is closed by a cone, the current remains 
finite. With the notation for the ends of each 
section as shown in figure (12), these conditions 
can be writ ten for the junctions as 

^pfPtt&i^) J p ( p n + |, aj B') 

(77a) 

j d(p/ib,B) — J e (pn+l,ajQ) 

(77b) 

By using equations (65) and (66), equal 
become 

ions (77) 

i ;; 1 1 It ft | j C rl [t Bfi - I + 

(78a) 

A n (l n jf B fl h n b~\~Jnb In -f- jR-ft + i . a + 

where 

Paa^+’-'i'Oifa 

n,ia 4— esc 2 <t> n 

a l 3 n 4-1 . a 

(78b) 

(79a) 

, Pw , cot 

4-cs c 2 <t> n 

(79b) 

Pna 00 1 <f) n 

“ 4— csc 2 0„ 

(79 c) 

2 p„„ cos <j>„ 
Jna 4— esc 2 4> n 

(79d) 

and a„„, b lltl] c„,„ and arc similarly defii 

icd. The 


two parts of equation (78) give, by addition and 
subtraction, 

\ . 1 C nb~\~jnb \ _ 1 Cn. + I,a”K/« + I,a fon«\ 

A n a nb + ~2 — =A n + 1 a a+lt a+- " 2 ~ l 8Ja / 
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B„b nt +**-^=B H+l b n+ i. a + C ” +L? - 2 ^ ±1 - (80b) 

In the application of these equations, it is im- 
portant to note that, in proceeding from one end 
of the frustum series to the other, if the opening 
angle of a cone with radius increasing with the 
distance along the body is considered positive, 
the opening angle of a cone with radius decreasing 
as z increases must be negative, and vice versa. 
This sign convention is intrinsically assumed by 
the equations. 

Equations (80) furnish 2m— 2 of the equations 
necessary for the determination of the 2m con- 
stants. Tn order to set up the remaining two 
equations, it is necessary to consider the following 
three cases, each of which must be treated 
separately: 

(1) Both ends closed: In this case, finiteness 
requires that 

7?! — 0 (81a) 

A m =0 (81b) 

Equation (81a) provides an initial value from 
which successive B n values can be calculated by 
the recurrence equation (80b). Likewise the A n 
values can be determined by equations (80a) and 
(81b). 

(2) One end open and one end closed: At the 
closed end, finiteness requires that 

/?! — 0 

and, at the open end, the requirement that the 
longitudinal current vanish gives equation (68b) 
applied to the mill segment. Tlius, the B n values 
can again be calculated by recurrence, after which 
equation (68b) gives A m , and the A n values can 
then be determined similarly. 

(3) Both ends open: The boundary conditions 
for this case are that the radial component of the 
current vanish at each end, so that equations (68a) 
and (68b) are applied to the first and last segments, 
respectively. Next A rn and B m are expressed by 
linear relations in A x and B u respectively, obtained 
by successive use of equations (80a) and (80b). 
Equations (68a) and (68b) now become a pair of 
simultaneous equations in A x and B i} so that these 
two can he found. The others then follow. 

The constants in the current equations having 
been determined, the flow is completely defined 
and the torque follows by summing the contribu- 


tions of all the sections, each being determined by 
equation (70). 

GENERAL THIN- WALL BODY OF REVOLUTION 

The next step is to apply the theory developed 
for series of frustums to the case of a continuous 
body of revolution (fig. 13) having as generatrix 
an arbitrary curve, say r=r(z), which lias a piece- 
wise continuous first derivative drjdz. Initially, 
however, it will he assumed that drjdz is continuous 
throughout the length of the body. 

The body is considered as made up of a largo 
number of frustums joined; therefore, equations 
(80) apply. Equation (80a) is now rewritten in 
the following form (for a reason that will he im- 
mediately apparent): 

A n i&n + 1,0 ’ A n d n t) 

= 2 K c nb—Cn+l.a)~\-(fnb—Jn+l.a)] (S2) 

If the body is considered to be continuously 
curved, 

1 — A | A ~ 

J 

t dd'n . ^ 

With these expressions, equation (82) becomes, in 
the limit, 

d[a(z)A(z)]=~ \dc(z)+dj{z)\ (S3) 

The subscripts arc no longer needed, as A, a, c, 
and / are continuous functions of z . Equation 
(83) may be integrated directly to give 

a(z)A(z) = ~ C - ^^ ) +a(0M(0) 

+ c(0)t M (S4 ) 




Fiofre 13 —Geometry of general body of revolution. 
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.1 (--)■=»-'( J> [ol0j.4{0) + c,0: ' ^ f " >! 

-»>] (85) 

Similarly, 

B{s)=b-\z) [i(Q)7?(0)+ c(0) +- /(0) 

_C(2)-/U)] (S6) 

The quantities a, b 7 c, and / are defined by equa- 
tions (79) as functions of p and 0. These are de- 
fined in terms of r(z) by the geometry of the prob- 
lem, as shown in figure (13), as 


/l A-( (lr \~ 2 

p V 1+ W 

0=tan _l (yjjzj 


(87) 


(SS) 


Other needed relations that follow from figure 13 
are 


2 , _ , (dr\ 2 

CSC 2 0 = 1+f ) 

fJL /dry 1 

cot 0=[ -7- 1 
\dzj 


The constants A (0) and /?(()) are determined in 
the same manner as that outlined for the joined 
frustums, in order to satisfy the condition that 
no current flows out the ends. The functions 
A(z) and B(z) thus evaluated arc substituted into 
equations (64), (65), and (66) to obtain the 
stream funetion and current components. The 
torque then follows as 


L—iraC 2 h 2 cok 


cos 2 0 sin 0 
Jp 4— esc 2 0 

a 


[p 3 +A(z)p C8C * +1 


+B(z)p- C8c * +1 ]dp (89) 


In general, this integral would have to be evaluated 
numerically. Any discontinuities in the slope 
arc accounted for merely by treating the discon- 
tinuity as a juncture in a series of frustums. 


NUMERICAL EXAMPLE 


micrometeoroid satellite S-55 in a symmetrical 
spin is calculated. For the calculations, the 
following numbers are used: 


h = 0.30 gauss 
co = 2 1 . 0 radian s/sec 
l— 96 cm 
r=24 cm 
7 = 0.05 cm 
c — 3 X 1 0 10 cm/sec 
a— 0.3 12X 10 6 (ohni-cm)” 1 
X — 90° 

2r l 

The fineness ratio is 2, for which 1 — tanh —=0.52. 

I 2 r 

Tt is necessary to express a in Gaussian units, 
that is, in (statohms-cm) -1 . The torque is then 
calculated by equation (19) as: 


L= 


— k7rX2.Sl X 10 17 X 


(Q.3) 2 X21.0 

9X10 20 


X24 3 X96X0. 05X0.52 


— 64.0k dync-cm 


RESUMfi 

A theoretical analysis has been made of the 
eddy currents induced by an applied magnetic 
field on the following spinning shapes: 

(1) Tliin-wall symmetrically spinning cylinder, 

(2) Thick- wall symmetrically spinning cylinder, 

(3) Thin-wall tumbling cylinder, 

(4) Thin-wall cone and cone frustum, 

(5) Joined thin-wall cone frustums, and 

(6) General thin-wall body of revolution. 

From the current expressions, the torques are 
calculated. The first two cases were solved by 
applying boundary conditions to the scalar potential 
of the electric field. The other cases were solved 
by means of a stream function. 

Figures that show the variation of torque with 
fineness ratio and thickness ratio are presented 
for thin- and thick-wall cylinders. From these, 
the degree of approximation in the thin-wall 
treatment can be ascertained. It was found that 
the average torque acting on a tumbling cylinder 
is one-half the torque acting on a symmetrically 
spinning cylinder, all oilier factors being equal. 


As an example of the application of the formulas 
of this paper, the magnetic torque acting on an 
aluminum cylinder such as the heat shield of the 
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